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We study nonlinear response in weakly coupled hot </> theory. We ob- 

. tain an expression for a quadratic shear viscous response coefficient using two 

\j3 i different formalisms: transport theory and response theory. The transport 

' bution function and expanding in the gradient of the local four dimensional 

■ velocity field. By doing a gradient expansion on the Boltzmann equation we 



obtain a hierarchy of equations for the coefficients of this expansion. 

To do the response theory calculation we use Zubrave's techniques in 
QhI nonequilibrium statistical mechanics to derive a generalized Kubo formula. 

Q^' Using this formula allows us to obtain the quadratic shear viscous response 

^ ■ from the three-point retarded green function of the viscous shear stress ten- 

sor. We use the closed time path formalism of real time finite temperature 
field theory to show that this three-point function can be calculated by writ- 
^ I ing it as an integral equation involving a four-point vertex. This four-point 

' vertex can in turn be obtained from an integral equation which represents the 

resummation of an infinite series of ladder and extended-ladder diagrams. 

The connection between transport theory and response theory is made 
when we show that the integral equation for this four-point vertex has exactly 
the same form as the equation obtained from the Boltzmann equation for the 
coefficient of the quadratic term of the gradient expansion of the distribution 
function. We conclude that calculating the quadratic shear viscous response 
using transport theory and keeping terms that are quadratic in the gradient of 
the velocity field in the expansion of the Boltzmann equation is equivalent to 
calculating the quadratic shear viscous response from response theory using 
the next-to-linear response Kubo formula, with a vertex given by an infinite 
resummation of ladder and extended-ladder diagrams. 
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I. INTRODUCTION 



Fluctuations occur in a system perturbed slightly away from equilibrium. The responses 
to these fluctuations are described by transport coefficients which characterize the dynam- 
ics of long wavelength, low frequency fluctuations in the medium The investigation 
of transport coefficients in high temperature gauge theories is important in cosmological 
applications such as electroweak baryogenesis ^ and in the context of heavy ion collisions 

i- 

There are two basic methods to calculate transport coefficients: transport theory and 
response theory [p|^T^. Using the transport theory method one starts from a local equilib- 
rium form for the distribution function and performs an expansion in the gradient of the 
four- velocity field. The coefficients of this expansion are determined from the classical Boltz- 
mann equation |]10|,|12[. In the response theory approach one divides the Hamiltonian into a 
bare piece and a perturbative piece that is linear in the gradient of the four-velocity field. 
One uses standard perturbation theory to obtain the Kubo formula for the viscosity in terms 
of retarded green functions 0. These green functions are then evaluated using equilibrium 
quantum field theory. As is typical in finite temperature field theory, it is not sufficient to 
calculate perturbatively in the coupling constant: there are certain infinite sets of diagrams 
that contribute at the same order in perturbation theory and have to be resummed |r3,H 



In this paper, we want to compare these two methods. The response theory approach 
allows us to calculate transport coefficients from first principles using the well understood 
methods of quantum field theory. On the other hand, the transport theory approach involves 
the use of the Boltzmann equation which is itself derived from some more fundamental theory 
using, among other things, the quasiparticle approximation. In this sense, the response the- 
ory approach is more fundamental than the transport theory method. However, the response 
theory approach can be quite difficult to implement, even for a high temperature weakly 
coupled scalar theory, because of the need to resum infinite sets of diagrams [jl3|,|l4|. These 



considerations motivate us to understand more precisely the connection between the more 
practical transport theory method, and the more fundamental response theory approach. 

Some progress has already been made in this direction. It has been shown that keeping 
only terms which are linear in the gradient expansion in the transport theory calculation is 
equivalent to using the linear response approximation to obtain the usual Kubo formula for 
the shear viscosity in terms of a retarded two-point function, and calculating that two-point 
function using standard equilibrium quantum field theory techniques to resum an infinite 
set of ladder diagrams |T0|,|15|,[T6| . This result is not surprising since it has been known for 
some time that ladder diagrams give large contributions to ra-point functions with ultra soft 
external lines 1113 . 



To date, calculations of transport coefficients have been limited to linear response. In 
some physical situations however nonlinear response can be important, especially for rela- 



tivistic gauge theories ||17|-[20|1. It is therefore of interest to study nonlinear response. In 
this paper we study nonlinear response to fluctuations in weakly coupled high temperature 
scalar 0^ theory using both the transport method and the response theory approach. We 
study the relationship between these two approaches at the level of quadratic response. 

This paper will be organized as follows. In section II we define shear viscosity and 
quadratic shear viscous response using a hydrodynamic expansion of the energy-momentum 
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tensor which includes up to quadratic terms in the gradient of the four-velocity. In section 
III we calculate viscosity using the transport theory method by performing a gradient ex- 
pansion on the Boltzmann equation and obtaining a hierarchy of equations. In section IV 
we derive the quadratic response Kubo formula by using Zubrave's techniques in nonequi- 
librium statistical mechanics. We obtain an expression which relates the quadratic shear 
viscous response coefficient to the three-point retarded green function of the viscous shear 
stress tensor. Starting from this generalized Kubo formula we calculate the quadratic shear 
viscous response using standard techniques of finite temperature quantum field theory. We 
show that the quadratic shear viscous response can be obtained as an integral over a four- 
point vertex. This four-point vertex satisfies an integral equation involving terms which are 
quadratic in a retarded three-point function, which itself satisfies a linear integral equation. 
These two integral equations represent the resummation of an infinite series of ladder and 
extended-ladder diagrams. We show that these integral equations, which represent ladder 
and extended-ladder resummations, have the same form as the first two equations in the 
hierarchy obtained from expanding the Boltzmann equation. We discuss our results and 
present our conclusions section V. 



II. VISCOSITY 

In a system that is out of equilibrium, the existence of gradients in thermodynamic 
parameters like the temperature and the four dimensional velocity field give rise to ther- 
modynamic forces. These thermodynamic forces lead to deviations from the equilibrium 
expectation value of the energy momentum tensor which are characterized by transport 
coefficients like the thermal conductivity and the shear and bulk viscosities. In order to 
separate these different physical processes we decompose the energy-momentum tensor as, 

T'^'^ = eu^v!" - pA'^^ + P^v!" + P^'u^ + tt'^'^ ; A^, = g^, - u^u,. (2.1) 

The quantities e, p, P^ and tt^j, have the physical meanings of internal energy density, 
pressure, heat current and viscous shear stress, respectively. The four vector u^{x) is the four 
dimensional four-velocity field which satisfies u^{x)u^{x) = 1. The expansion coefficients 
are given by 

p=-^A«^T"^ 
P^ = A^^u^T-f' (2.2) 
TT^, = (A^„A,^ - ^A,^A„^)T-^. 

The viscosity is obtained from the expectation value of the viscous shear stress part of 
the energy momentum tensor. We expand in gradients of the four-velocity field and write, 

d{n,^) = rj^'^H,^ + r^^'^Hl^ + • • • (2.3) 

2 
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where rj^^^ and r]^"^^ are the coefficients of the terms that are hnear and quadratic respectively 
in the gradient of the four-velocity. The first coefficient is the usual shear viscosity. The 
second has has not been widely discussed in the literature - we shall call it the quadratic 
shear viscous response. 

Throughout this paper we work with 0^ theory. The Lagrangian for this theory is 

^ = ^[(5.0r-^V]-^0' (2.4) 
and the coupling constant is assumed to be small: A ^ 1. 

III. VISCOSITY FROM TRANSPORT THEORY 

Kinetic theory and the Boltzmann equation can be used to calculate transport prop- 
erties of dilute many-body systems. One assumes that, except during brief collisions, the 
system can be considered as being composed of classical particles with well defined position, 
energy and momentum. This picture is valid when the mean free path is large compared 
with the Compton wavelength of the particles. At high temperature the typical mean free 
path of thermal excitations is (9(1/A^T) and is always larger than the typical Compton 
wavelength of effective thermal oscillations which is 0{1/VXT) [|T^]. We introduce a phase 
space distribution function f{x,k) which describes the evolution of the phase space proba- 
bility density for the fundamental particles comprising a fluid. In this expression and in the 
following equations the underlined momenta are on shell, since we are describing a system 
of particles. The form for f{x,k) in local equilibrium is, 

/^°^ = ,/.(.K(I).^_l ^=^^; N,:=l + 2n,. (3.1) 

We study the Boltzmann equation in the hydrodynamic regime where we consider times 
which are long compared to the mean free time and describe the relaxation of the system 
in terms of long wavelength fluctuations in locally conserved quantities. For a simple fluid 
without any additional conserved charges, the only locally conserved quantities are energy 
and momentum. To solve the Boltzmann equation in this near equilibrium hydrodynamic 
regime, we expand the distribution function around the local equilibrium form using a 
gradient expansion. We go to a local rest frame in which we can write u{x) = 0. Note that 
this does not imply that gradients of the form diUj must be zero. In the local rest frame 
( p.3|) becomes, 

S{7c,,) = -v^''>H,,+r^^''>Hj:' + --- (3.2) 
2 

Hij = diUj + djUi - -6ij{d ■ u) 

Hjj^ := Hi^H^- — -6ij HimHim 

In all of the following expressions we keep only linear terms that contain one power of Hij 
and quadratic terms that contain two powers of Hij, since these are the only terms that 
contribute to the viscosity coefficients we are trying to calculate. 
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We write, 

/ = /(°) + f^'^ + f^'^ + ■■■ (3.3) 

with, 

/W-i^S'^/W; (3.4) 

Using (|3.1| ) we obtain, 

k^d'^fW = -f3nk{l + n,) \h,{k)H,j (3.5) 
/(I) := -nk{l + nk)(t)k ; 0fc = P^Bij{k)Hij 

:= nfc(l + nk)Nueu ; := /3^]Q,Uk)H,jHi 



where we define 



and write, 



4(^) = (^.%-^^..); h,{i^ = i^k{k) (3.6) 



5(A;)i, = iim{k)B{k) ■ Ci.Uk) = ii,{k)iUk)C{k). (3.7) 
The viscous shear stress part of the energy momentum tensor is given by 

d^k 



Using the expansion (^]3|) we get. 

The lowest order term is identically zero. We obtain the linear and quadratic contributions 
by substituting in ( |3.5| ) and (p.7|). We use Ha := and the following results which are 
obtained from rotational invariance: 

kik,B{k) ^ ^6^,eB{k) 

kikjkik^B(J^ ^ Y^('^«i'^'™- ~^ ^u^jm 

+ Si^Sji)k^B{k) 

kikjkikrnkakbC{k) — [ 6ab{SlmSij + Slj6mi + 5ii5mj) + 5a/(5fem% + 6biSmj + hj^mi) 

+ 5aj(5w5mj + ^bm^ii + ^U^ml) ] k'^C{k) (3.10) 
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We obtain, 



Comparing with ( p. 21) we have, 

^^'^ = ^ / ..1?. Ml + n,)eB{k) (3.11) 



15 7 (27r)32cjfc 

Thus we have shown that the shear viscosity and the quadratic shear viscous response can 
be obtained from the functions B{k) and C{k) respectively. These two functions are the 
coefficients of the hnear and quadratic terms in the gradient expansion of the distribution 
function. In the next section we will show that these functions can be obtained from the 
first two equations in the hierarchy of equations obtained from the gradient expansion of 
the Boltzmann equation. 



A. Expansion of the Boltzmann Equation 

The Boltzmann equation describes the evolution of the distribution function f{x,k) and 
can be used to obtain integral equations for the functions B{k) and C{k) defined in ( p.5|) 



and (|3.71 ). The Boltzmann equation has the form: 

k^d^f{x,k)=C[f] (3.13) 

where C[f] is the collision term: 

C[f] = 1 [ rfr«[/l/2(l + /3)(1 + fk) - (1 + /l)(l + /2)/3/fc] (3.14) 
2 J123 

with fi := f{x,p.), fk '■= f{x,k). The symbol dri2^3k represents the differential transition 
rate for particles of momentum Pi and P2 to scatter into momenta P3 and K and is given 
by 

dTu^,k := ^imP3'P2'Eil'nii^^^^^(27r)^5(p^ + -P3 - Ai) (3.15) 



where T is the multiparticle scattering amplitude. Using the expansion (|3.3|) and ( p.5|) 
produces a hierarchy of equations. 

1. First Order Boltzmann Equation 
The first order equation is: 

!fd,f{x,k)=C[f^'^-J^'^] (3.16) 
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where we keep terms linear in f^^^ on the right hand side. The left hand side gives, 

]fd^f{x,k) = -nk{l + nk)(3h,,{k)Hij (3.17) 

The right hand side is 

C[f^'^] = \ I ^^ri2^3fe(l + ^i)(l + n2)n3nfc(0fc + 03-0i-02) (3.18) 

Z J 123 

Using the definition of given in (|3.5| ) and comparing the coefficients of Hij on both sides 
of ( p.l6| ) we obtain, 



h^k) = \l dV,,^,,^-^±^^^^±^[B,,{p;) + B,,{pJ - B,,{k) - B,,{pJ ] (3.19) 

This result is an inhomogeneous linear integral equation which can be solved self-consistently 
to obtain the function Bij{k). 

2. Second Order Boltzmann Equation 

The second order contribution to ( |3.13| ) is, 

ra^/('n^,fc)=C[/(°);/(^);/(^)] (3.20) 

where we keep terms linear in f'^'^^ and quadratic in f^^^ on the right hand side. Using ( p.5|) 
the left hand side becomes, 

fcM^^/d) = nu{l + nu)Nul3^hi,{k)Bi^{k)H,,Hi,r. 
and the right hand side gives, 

^(2)^^(0). ^(1). ^(2)] ^ 1 r dT,2^^,{l + n,){l+n,)n^nk 

I J 123 

[(iVi^^i + N^e^ - Ns9s - N,9k) 

+ ^{{Ni + N2)M2 - {Nk + N3)^3(Pk + (A^3 - A^l)0103 

+ {Nk - N,)(l),(t>k + (iVa - iV2)0302 + (iVfc - iV2)0fc02)] (3.21) 

Using the definitions of (p and 6 given in ( |3.5|) and comparing the coefficients of HijHim on 
both sides we obtain, 

nk(l+nk)NkIij{k)Bim{k) = ^ [ dT i2^3k{^ + ni)(l + n2)n3nk 

2 J123 

{ [NidjiUp,) + N2Cijim{p^) - NkCiji^it) - NsCiji^p^)] (3.22) 

+ i[((iVi + N2)B,,{p^)BiUp2) - {Nk + N,)B,,{pjBUk) + {N^ - N,)B,,{p^)BiUp,) 
+ {Nk - m)B,,{p^)Bimik) + (iVg - N2)B,,{p^)BiUp^) + {Nk - N2)B,^{k)BiUp^)]} 

This integral equation can be solved self consistently for the quantity Cijim{k) using the 
result for Bij{k) from (|3.19| ). 



IV. VISCOSITY FROM FIELD THEORY 



The Kubo formulae allow us to use quantum field theory to calculate nonequilibrium 
transport coefficients. The results should be the same as those obtained in the previous 
section using transport theory. To do calculations in a system that is out of equilibrium, the 
Hamiltonian is separated into an equilibrium piece Hq and a nonequilibrium piece Hext which 
depends on the gradients of the thermodynamic parameters: the four-velocity field and the 
inverse temperature field. For systems close to equilibrium the nonequilibrium piece of the 
Hamiltonian can be treated as a perturbation and the deviations of physical quantities from 
their equilibrium values can be calculated perturbatively. The linear response calculation 
includes only the first order contribution to this expansion and gives transport coefficients 
that can be expressed as integrals of retarded two-point Green functions over space and 
time. One of the results we obtain in this paper is that the quadratic shear viscous response 
coefficient can be written in terms of a retarded three-point function. Throughout this 
section we use capital letters to denote four-vectors and small letters for three-vectors. We 
also define J d^p/ {27r)^ := f dP. 

To calculate the expectation value of an operator we take the trace over the density 
matrix. We follow the presentation of We work with the density matrix in the Heisenberg 
representation which satisfies, 



and can be written as. 



Tre-^+-s 



(4.2) 



where 



A = I d^xF^Tr 



t 



B=fd^xf dt'e'^''-'^T^^{x,t')df'F''{x,t') (4.3) 

J J —oo 

with F^ = f3u^ and e to be taken to zero at the end. In this expression A is the equilibrium 
part of the Hamiltonian and B is a. perturbative contribution that is linear in the gradient 
of the four-velocity field. Note that in the local rest frame = (1,0,0,0) and [A]^^-^ = 
/ d^x PTqq. Using the identity for the exponential function of two operators 

_ g/3aM ^ dXeHe-^^^+''^] (4.4) 
Jo 

we can expand 

^-A+B ^ g-AM ^ dXe^A-^^-XA 

Jo 

+ f dX f dre^^Be-^^e^^Be-^^ + C(5=^)] (4.5) 
Jo Jo 
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Thus the density matrix can be written as 

p = po [ 1 + /' d\{B{\) - {B{\))) + /' d\ r dTiB{\)B{r) - {B{\)B{r))) 
Jo Jo Jo 

- fd\ ['dX'{{B{X))B{X')-{B{X)){B{X'))] + 0{B') (4.6) 
JO Jo 



where 

is the local equilibrium density matrix. 

A. Viscosity as an Expansion in Green Functions of Composite Operators 

1. Linear Response 

Using the first three terms of ( ^.6| ) produces the linear response approximation for the 
deviation of the expectation value of the viscous shear stress part of the energy momentum 
tensor from the equilibrium value: 

'rfV r dt'e'^''-'\7iap{x',t'),7T^^{x,t))d"F^{x',t") (4.8) 



where the correlation function {7io,f3{x',t'),7if^,^{x,t)) is defined as 

pJo 



1 fl^ 

7r^^{x, t), Tla/six', t')) = - dT{{'Kap{x', t' + Zr)7r^^(x, t) - {Hafsix', t' + iT))QTXf,^{x, t))o 



(4.9) 



and (■ ■ ■)() := Tr(po ■ ■ ■)• From now on we drop the subscript on the correlation functions. 
We take the initial time to to minus infinity and we assume the system is in equilibrium at 
t = to and that the external forces are switched on adiabatically. We assume that changes 
in the thermodynamic forces are small enough over the correlation length of the correlation 
functions that factors = (3u^ can be taken out of the integral. Using ( p.2| ) and rotational 
invariance we have, 

1 2 

^^^lv'^a|3 = — (A^oAj,^ + A^/^A^^ - -^fj.u^af3)'n'aT tt'^'^ (4.10) 

Using this result we obtain, 



10 

13 

{{{n{x', t + zr)7r(x, t)) - {{n{x\ t + ir))7r(x, t))) (4.11) 
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where we have written tTo-tTt'^'^ := vrvr. If we assume that correlations vanish as t' —oo 
this expression can be rewritten as, 

r r dt"-^{{{n{x',t" + ZT)7i{x,t)) - {{7r{x',t" + ^r))7r{x,t))) (4.12) 
Jo J —oo dz 

Using -^f{t" + ^t) = —i^f{t" + ir) we can perform the integration over r. We obtain, 
6{7r^.y (4.13) 
= -^H^, f d?x' f /* dt"((7r(x',t" + ^/5)7r(x,t)) - (7r(x', t")7r(x, t))) . 

lU J J —oo J —oo 

Using the KMS condition 

(7r(x', t" + tl3)7T{x, t))= Tr[e-^^7r(x', t" + i(3)tt{x, t)] 

= Tr[e-'^^7r(x, t)7r(x', t")] = (7r(x, t)7r(x', t")) (4-14) 



we obtain 



lU J J —oo J —oo 

where 



Dr{x, t- X, t") = -ie{t - t") [7r(x, t), 7r(x', t")] (4-16) 
We extract the shear viscosity using the definition ( p.3[ ). We obtain, 

^(1) = L( f dt'e<''-'^ f dt"DniO;x',f) (4.17) 

lU J J —oo J —oo 

We can rewrite this result in a more useful form. We rewrite the integrals inserting theta 
functions: dt' = 0{—t') dt' and use the integral representation for the theta function: 

0(,)^I^J:^. (4.18) 

Inserting D^j(0; a;', t") = J dQ e'^^' Dn{Q) with X' = {x',t") we obtain. 

Since it is the imaginary part of the two-point function that is odd in go we have, 

V^'^ = ^-^lm[]imDnmUo=o (4.20) 
10 ago 9-^0 

This is the well known Kubo formula which expresses the shear viscosity in terms of a 
retarded two-point green function that can be calculated using equilibrium quantum field 
theory. 
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2. Quadratic Response 

Now we consider corrections to the linear response approximation. We calculate the 



quadratic shear viscous response from the terms in (|4.6|) that are quadratic in the interaction. 
We show that the result can be written as a retarded three-point correlator. We obtain, 

rt . . . rt 



J J J ~oo J —oo 

(7r^,(x, t), n^pix', t'), TTp^ix", t"))d'^F^{x\ t')d''F''{x\ t") . (4.21) 



We use (2.2) and rotational invariance to write, 



■■= ^M^'^ n) (4.22) 

Writing out the correlation function we obtain, 

6Mx,t)y = ^Hl^ I rfV / d^x" f dt'e<''-'^ f rft"e^(*"-*) 

OO J J J — oo J —oo 

[ dr r d\{7i{x',t' + iT)7r{x",t" + i\)7i{x,t)) 
Jo Jo 

— I dr dr' {{^{x^t' + iT))n{x" ,t" + iT')n{x,t)) 
Jo Jo 

dr I d\{{-K{x',t' + iT)Tc{x",t" + iX))Tc{x,t)) 
10 Jo 

+ dr dT'{{n{x',t' + iT)){7r{x",t" + iT'))7r{x,t))] (4.23) 
JO Jo 

Assuming that correlations vanish when the time approaches minus infinity we can rewrite 
O) as 

5{n^,{x,t))'' = —Hl^ I d'x' I d^x" f rft'e^(*'-*) /* dt"e<'"-''^ f ds' f " ds" 

35 J J J —oo J —oo J —oo J —oo 

-T-i-rrl ['^ dr I dX{n(x', s' + iT)n(x", s" + zA)7r(x, t)) 
ds' ds' Jo Jo 

dr / dT'{{-n{x' ,s' + iT))-n{x" ,s" + iT')Ti{x,t)) 
lo Jo 

— / dr I dX{{iT{x' , s' + iT)-K{x" , s" + iX))TT{x,t)) 
Jo Jo 

+ f dr dT'{{TT{x',s' + tT)){n{x",s" + tT'))n{x,t))] (4.24) 
Jo Jo 

Carrying out the integration over A, r and r' and using the fact that we have symmetry 
under interchange of {x',t') and {x",t") we obtain, 

5(7r^,(x,t))^ = f d^^' f ^3^// /•* dt'e^^''-'^ f rft"e^(*"-*) /*' ds' f" ds" 

i\j J J J —OO J —OO J —oo J —oo 

^iMx,t),nix',s%nix",s")] + [[7r(x, t), 7r(x", .")], vr(x', .')]) (4.25) 
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We define 

t; x\ s'- x", s") = ^{[[n{x, t), n{x', s% 7i{x", s")] + [[7r(a;, t), 7r{x", s'%rr{x', s')]) (4.26) 

and write 

6{n^,{x,t)y 

TO J J J — OG J — OO J — OO J — CO 

Using (|2.3| ) we extract, 
3 
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/ d^x' r dt'e''' r dt"e''" f ds' f ds"G{<d; x\ s'; x" , s") (4.27) 

J J—OO J — OO J—OO J—OO 



As we did previously, we can rewrite this result in a neater form. We rewrite the integrals 
inserting theta functions: /^o^rft' = 6{—t') J^^dt' and use the integral representation for 
the theta function ( [4.18| ). We use X' = (x', s') and X" = {x", s") and rewrite, 

G(0; x', s'; x\ s") = J dQe''''^ J dg's'"""^' G{-Q - Q' , Q, Q') . (4.28) 

We obtain. 



Using the techniques of to write the three-point function in the spectral representation, 
it is tedious but straightforward to show that the only contribution to the result above 
comes from the real part of the three-point function which is retarded with respect to the 
first leg. In coordinate space this retarded three-point function is written: Gjii{x,y,z) = 
e{U - ty)e{ty - Q[[n{x), 7r(y)], tt{z)] + e{U - Qe{t, - ty)[[Tr{x), 7r(z)], 7r(y)]. We obtain. 

This is an interesting new result. We have obtained a type of nonlinear Kubo formula that 
allows us to obtain the quadratic shear viscous response from a retarded three-point function 
using equihbrium quantum field theory. 



B. Diagrammatic Expansion 

We obtain a diagrammatic expansion for the viscosity coefficients given in (|4.20|) and 
( [4.30| ). We use the closed time path formulation of finite temperature field theory, and 
work in the Keldysh representation. Several reviews of this technique are available in the 
literature p2|-p6|. The closed time path integration contour involves two branches, one 
running from minus infinity to positive infinity just above the real axis, and one running 
back from positive infinity to negative infinity just below the real axis. All fields can take 
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values on either branch of the contour and thus there is a doubhng of the number of degrees 
of freedom. It is straightforward to show that this doubhng of degrees of freedom is necessary 
to obtain finite green functions. We discuss below the structure of correlation functions of 
field operators. 

The two-point function or the propagator can be written as a 2 x 2 matrix of the form 

where Dn is the propagator for fields moving along Ci, D12 is the propagator for fields 
moving from Ci to C2, etc. The four components are given by 

Dn(X-F) = -^(T(0(X)0(r))) 
D^2{X-Y) = -l{<P{Y)<|){X)) 
D2i(X-F) = -^(0(X)0(F)) 

D22{X -Y) = -i{f{<P{X)<P{Y))) (4.32) 

where T is the usual time ordering operator and T is the anti-chronological time ordering 
operator. Physical functions are obtained by taking appropriate combinations of the com- 
ponents of the propagator matrix. It is straightforward to show that the usual retarded and 
advanced propagators: Dji = —i9{xo — yo)[(j){X), (j){Y)] and Da = —iOiyo — Xq)[4>{X), 0(^)] 
are given by the combinations, 

Dr = Dn - Du 

Da = Dn - D21 . (4.33) 

The IPI part of the two-point function, or the polarization insertion, is obtained by trun- 
cating legs. The retarded and advanced parts are given by, 

= Un + . (4.34) 

The situation is similar for higher n-point functions. For example, the three-point func- 
tion which is retarded with respect to the first leg is given by 

Tm = Tin + rn2 + ri2i + ri22 ■ (4.35) 

The other three-point vertices that we will need are: 

= Till + rii2 + r2ii + r2i2 
r_R3 = Till + ri2i + r2ii + r22i 
Tf = Till + ri2i + r2i2 + r222 

The four-point function which is retarded with respect to the first leg is given by, 

Mri = Mini + M1112 + M1121 + M1211 + M1122 + M1212 + M1221 + M1222 • (4.36) 
The other four-point vertices that we will need are: 
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Mm = Mini + M1121 + M1211 + M2111 + M1221 + M2121 + M2211 + M2221 

Mp = Mini + M1121 + M1211 + M1221 + M2112 + M2212 + M2122 + M2222 • (4.37) 



There are two relations that we can use to simphfy expressions involving these vertices. The 
first is a consequence of the fact that real time green functions are related to each other 
through the KMS conditions. These identities are a consequence of symmetries generated 



through cyclic permutations of field operators. The identity that we will need is 21 



+ NiTrs + iVgr^ = (iVi + N3)r^^ (4.38) 

where each vertex carries arguments r{Pi, P2, P3) and we write N{Pi) := A'^i etc. The 
second simplification occurs because of the fact that the equations we obtain contain a 
specific combination of four-point vertices which we define as, 

Mf = Mf + NiMr^ + N^Mri (4.39) 

where each vertex carries arguments M(Pi, P2; Pzi -^4)- It has been shown that this particular 



combination of four-point vertices has special properties in a number of other contexts [27 

In addition, bare vertices that carry a Keldysh index '2' have an extra factor of minus 
one. This factor is accounted for by inserting a factor of r for each vertex where r is the 
two component vector given by. 




We want to obtain a perturbative expansion for the correlation functions of composite 
operators Dr{X, Y) and Gri{X, Y, Z) which appear in (|4.2(]| ) and ([4.30|) . We obtain expres- 



sions for these correlation functions in terms of the vertices Fjj and Mijim which are defined 
as the vertices obtained by truncating external legs from the following connected vertices: 

rg = (T,7r,,.(x)0(r)0(z)) 

tC 



Mtjim = (r,7r,,(X)7r,„(r)0(Z)0(W^)) (4.40) 



where 



7T,,{X) = d,<f){X)d,<j){X) - ^%(9„0(X))(9™0(X)) (4.41) 

and Tc is the operator that time orders along the closed time path contour. These definitions 
allow us to write the two- and three-point correlation functions as integrals of the form. 
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(a) 




(b) 

Fig. [1]: (a) Two-point function for shear viscosity from linear response; (b) 
Three-point function for quadratic shear viscous response. The dashed external 
line represents the composite operator nij. The square box is the four-point 
function M and the round blob is the three-point vertex T. 



DabiQ) = 2ij dPTlUP, Q, -P - Q)iD,,{P)iD,,{P + Q)I,,{p)n (4.42) 
GaU-Q - Q', Q, Q') = ^j dK Mii^,T,tD,,{K)tD,a{K + Q + Q')Iu{k) 

where the indices {a, 6, c, rf, e} are Keldysh indices and take values {1,2}. These expressions 
are shown diagrammatically in Fig. [1] . We perform the sum over Keldysh indices using the 
Mathematica program described in [^]. We obtain, 

Dr{Q) = I J dK{Nk+, - N,)r%,{K, Q, -K - Q)Da{K)Dr{K + Q)I,,{k) (4.43) 
Gri{-Q-Q',Q,Q') (4.44) 
= -Aj dK{MF)ikkj{K, Q, g', -K-Q- Q')Dr{K)Da{K + Q + Q')I,,{k) 

To rewrite these expressions in a simpler form we use rotational invariance to write, 

Mijim := hjiimM; F^, := 4f (4.45) 

Using an obvious notation we write the pairs of propagators Da{P)Dr{P + Q) and 
Dr{K)Da{K + Q + Q') as ttpTp^q and ri^ak^q^qi . When {q'g , q'^} — > the dominant contribu- 
tion to the integral from these pairs of propagators is produced by what is called the pinch 
effect: the contour is "pinched" between the poles of the two propagators which gives rise to 
a factor in the denominator that is proportional to the imaginary part of the propagators. 
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We regulate the pinching singularity with the imaginary part of the hard thermal loop self 

(4.46) 



energy and obtain [|15 



Pk 

'2ImSL 



Pk = i{rk - ak) 



where E is the retarded part of the hard thermal loop self energy. 

Now we expand in go and ^q. In ( 4.43| ) we keep terms proportional to go since these terms 
are the only ones that contribute to ( |4.20| ); in ( [4.44| ) we keep terms proportional to gogp since 
these terms are the only ones that contribute to (|4.3(]| ). In each term there are products 
of thermal factors of the form A^^; — N^+q and A^^ — Nx+q' where x is some combination of 
{fc, p, r}. The expansion of these thermal factors is straightforward: 

Nx - Nx+q = 2go/5n,(l + n,) + ■ ■ ■ (4.47) 

Consider the behaviour of the vertices when {go, go} 0. Using ( [4.43|) and ( |4.46| ) in 
( [4.2(J| ) it is easy to see that only the real part of Tii2 contributes. Similarly, using ( |4.44D 
and ( |4.46| ) in ( [4. 301) it is clear that we need only the real part of Mp. First we consider 
( [4.43| ). Because of the explicit factor of {Nk+q — N^) we can set Q to zero in the vertex. 
The expansion of ( [4.44| ) in Q and Q' is more difficult. One can show that the terms in the 
expansion that contain gradients acting on the vertices Mp, M^i and M^4 do not contribute 
to the result and that we can make the replacement: 



ReMpiK, Q, g', -K-Q-Q')^ {N^ 



Nk+q+q')ReMm{K, 0, 0, -K) (4.48) 

= T{K) and 



From now on, to simplify the notation, we define T{K,0,—K) : = 
M{K,0,0,-K) := M{K). 

Using these results to simplify ( |4.43| ) and ( [4.44| ) and substituting into ( ^4.20] ) and ( [4.3U[ ) 
we obtain, 
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J dKk'^ Pknk{l + TLk) 



RerR2(K) 



(2) 



2^ 
105 



dKk PkUkil + nk)Nk 



ImSk 

ReMRi(K) 



ImSb 



(4.49) 
(4.50) 



Comparing with ( p.ll|) and ( |3.12| ) we see that the results are identical if we identify 

RerR2(k) 



B{k) 
C{k) 



ImSk 
ReMRi(k) 

ImSk 



(4.51) 
(4.52) 



with the momentum K on the shifted mass shell: 6{K'^ 



m 



th) 



where miu = rn? + ReT^K- 



1. The Ladder Resummation 

It has been known for some time that the set of diagrams which give the dominate 
contributions to the vertex Tij{P,Q, —P — Q) are the ladder diagrams. These diagrams 
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contribute to the viscosity to the same order in perturbation theory as the bare one loop 
graph and thus need to be resummed This effect occurs for the following reason. It 



appears that the ladder graphs are suppressed relative to the one loop graph by extra powers 
of the coupling which come from the extra vertex factors that one obtains when one adds 
rungs (vertical lines). However, these extra factors of the coupling are compensated for by a 
kinematical factor. This factor arises through the pinch effect which is described above. The 
addition of an additional rung in a ladder graph always produces an extra pair of propagators 
of the form akTk+q or r^afc+g. Products of this form contribute a factor which produces an 
enhancement. This factor occurs when the contour is "pinched" between the poles of the 
two propagators, which gives rise to a contribution in the denominator that is proportional 
to the imaginary part of the inverse propagators. 

In order to include ladder diagrams we obtain the vertex Tij{P, Q, —P—Q) as the solution 
to the integral equation shown in Fig. [2]. 





+ 1/2- 




Fig. [2]: Integral equation for the ladder resummation. The blob represents the 

vertex F. 



= I dPdRT',T^T^IiMDcaiR)DaciK + R- P)D,,{P + Q)Dab{P) (4.53) 

+^ j dPdRP,Z{P, Q, -P - Q)DUP)D,,{P + Q)D,,{R + K- P)T^D,^{R)r^ 



We perform the sums over the Keldysh indices using the Mathematica program in and 
simplify the result by taking Q to zero, keeping only the pinching terms, and using ( [4.46 ). 



The vertex function which includes the tree vertex is obtained by shifting: Fj^ = lij + ^F^-. 
(The factor of 1/2 is a symmetry factor). We obtain, 

-plm I 



1 

— C 

3 



k„^ki - -6mik'^ - — / dPdRdP' 



A2 

T 



'rlm(p\ 



(l + r^p)(l+r^p,)r^r/(l + ^fc) (4.54) 



Note that this integral equation is decoupled: the only three-point vertex that appears is 
F^. To simplify this expression further we use (|4.45| ), (|4.51| ) and the fact that Tfi2{P) is 
pure real, and symmetrize the integral on the right hand side over the integration variables 
{P, P', R}. We multiply and divide the left hand side by ImE^ and replace this expression 



in the numerator by the HTL result 10 
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J dP dR dP' {2n)^6\P + P' -R-K) 



^(^' ^)'™ = T^TTT^T ^ / dP' {2'n)%P + P'-R- K)p^p,p^, (4.56) 

12 1 + Uh) J 



ImSi. = — — ( 

12Vl + nfc. 

PpPrPp'{l + np)(l + np')nr (4.55) 
Rearranging we obtain |ll5| , 

A' 

12(1 +nfc) 



where we have used ( p.7| ), ( [4.45| ) and ( [4.51| ). When the delta functions are used to do 



the frequency integrals, this equation has exactly the same form as the equation obtained 
from the linearized Boltzmann equation ( ^.19| ) with a shifted mass shell describing effective 
thermal excitations. Comparing ( p.llD and (p.l9|) with ( [4.49|) and (|4.56|) we conclude that 



calculating shear viscosity using effective transport theory by keeping only first order terms in 
the gradient expansion is equivalent to using the Kubo formula obtained from linear response 
theory, with a three-point vertex obtained by resumming ladder graphs. The contributions 
to the viscosity are shown in Fig. [3]. 







Fig. [3] : Some of the ladder diagrams that contribute to shear viscosity. 



2. The Extended- Ladder Resummation 

In this section we consider an integral equation for the vertex Mijim which resums an 
infinite set of graphs that includes ladder graphs and some other contributions which we 
will call extended ladder graphs. We will show that this integral equation has exactly the 
same form as the integral equation ( p.22| ) obtained by expanding the Boltzmann equation 
to second order. We consider the integral equation shown in Fig. [4]. 




(») (l-) (c) 

Fig. [4]: Integral equation for an extended-ladder resummation. 
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Following (^4.39|) and ( |4.44|) we calculate contributions to Mp-, Mm, M^^ for each dia- 
gram. We introduce some notation to simplify the equations. For four-point vertices we list 
the first three momenta only and for three-point vertices we list the first two only. In both 
cases the last momentum is the one that is determined by the conservation of energy. For 
example: M{K, Q, Q\ -K-Q- Q') := M{K, Q, Q') and T{K, Q, -K - Q) := T{K, Q). In 
addition we write (as before) M(/s:, 0, 0, -fsT) := M{K) and T{K,{],-K) := T{K). From 
Fig. [4a] we have, 

^Mf^f-{K, Q, Q') = 2i[rk+rL{K, Q)T%{K + Q, Q') 

+ au+,V%,{K, Q)V%{K + Q') + h+,Tl,{K, Q)T'S:,{K + Q, Q')] 
{MUt^riK, Q, Q') = 2tn^,r%,{K, Q)Tt,{K + Q, Q') 

{M'itriK, Q, Q') = 2za,^,r%,{K, Q)T%{K + Q') (4.57) 

We can rewrite this result by substituting in the expanded form of the vertex Vij — lij + ^^ij 
with Tij obtained from the integral equation that corresponds to Fig. [2]. We use, 

= r^{k) -'-^JdPdR r,a,^,{N' - iV,,)a,r' [TUP, Q) - N,+,Vt,{P, Q) + N,V%{P, Q)] 
= P^ik) -'-^JdPdR r,a,+,iNr - iV')r,a' [TUP, Q) ~ K+A^P, Q) + Kn^^P, Q)] 

= I dPdR rpap+,(l - N'N,)prp' [V%{P, Q) - N,^,Tl,{P, Q) + N,V%{P, Q)] 

It is easy to show that the contribution from Fig. [4a] is equivalent to the contributions from 
the four diagrams in Fig. [5]. 



1/2- 



(a) (b) 




Fig. [5]: The four-point vertices that correspond to the diagram in to Fig. [4a] 



From Figs. [5a,5b] we obtain, 

{Mf^f^iK, Q, Q') = th,{k)[au+,T%{K, Q) + h^,T'^,{K, Q)] 
{M'l^'f'^iK, Q, Q') = in^,k,{k)T'^,{K, Q) 
{Mf,^f\K, Q, Q') = iak+M^)^m{K, Q) 



(4.58) 
(4.59) 
(4.60) 
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(Aff^f\K, Q, Q') = ih^{k)[Tk+qT%{K + g, Q') + h+qr%{K, Q)] 
{M'i'^f'iK, Q, Q') = ir,+qIij{k)r^SliK + Q, Q') 
{M'it)'\K, Q, Q') = ia^^qk^{k)T%{K + Q, Q') 



{Mii'"')'%K,Q,Q') 

= J dP dR ap+q+g>rp 

[^'f2{p + Q, Q') - A^p+,+,'r^i(P + Q, Q') + Np+qV%{p + g, g')] 

• [rF2(-^> Q){frap' + rrfp')ak+q 

+ ^Rli^, Q){fk+q{frap' + Trfp') + Vk+qUp'fr + Op'^r- + Va^))] 



[r^l(P + Q, QO - Np^q+q,T'^{{P + Q, Q') + Np+qT%{P + Q, g')] 



[r^2(^ + QO - Np+q+q,T'^^{p + g, g') + Np+qT'^^{p + g, g')] 

• rfl3(-ftr, Q)ak+q{fp'rr + ap'fr). 

{Mf'^f\K,Q,Q') 

= c/Prfi?a,+,r,[r'™(p,g)-iv,+,r^™(p,g) + iVpr^™(p,g)] 
[r^2(^ + Q,Q'K+aK/' + /y) 

(ii4^)5^(ir,g,g') 

= ciPciPrpap+,[r^-(^,i3)-^p+rt(p,g) + 7Vpr^(p,g)] 

• ri(K + g,g')rfe+,(a,/' + /y) 

(Mi^"^)^'^(x,g,g') 

= rfPrfPw,[r'™(P,g)-iVp+rt(p,g) + Ar^r^l(p,g)] 
■ r)^3(ir + g,g')afc+,(a7, + /v,) 



Prom Figs. [5c, 5d] we obtain. 




dPdR 



The contribution from Fig. [4b] is 



20 



^Mf^Y\K, Q, Q') = -y / dPdRap+,+,,rpMf^{P, Q, QOl^rV + a,v + Mp'] (4.70) 

{Miil"'y\K,Q,Q') = -y / dPdRa,+,+,,r,Mf^{P,Q,Q')[arU> + /.v] (4.71) 

(Mi^f )^''(ir,Q,Q') = -y / dPdRap+,+,,rpMf'^{P,Q,Q')[rrU, + f^a,,] (4.72) 
where P' = R + K — P. Fig. [4c] gives, 

(M^^'™)^'^(if,Q,g') = ^J dPdR{Np - N,^,){Nr-,^ - K,) (4.73) 
{M'l^Y\K,Q,Q') = ^J dPdR{N, - N,^,){Nr.q^ - Nr)r^{P,Q) (4.74) 

^ R2 (-^ Q iQ ^^r'^p^r+q^r^r—q' 

{Milt)''iK,Q,Q') = dPdR{N, - iV,+,)(iV,_,, - N^T^{P,Q) (4.75) 

^ R2 Q 1 Q ^Q'p'^pOjr+qQ'r'^r-q' 

To combine these expressions we use ( ^.38| ) and ( |4.39| ). We keep only the pinching 
contributions and use ( [4.46| ) to regulate. As discussed previously, we expand in go and 
and keep the term proportional to q'oQ'o' since that is the only term that will contribute to 
the quadratic shear viscous response coefficient. As before, one can show that the terms 
in the expansion that contain gradients acting on the vertices Mp-, Mm and do not 
contribute to the result and that we can use ( |4.48|) . We simplify further by using ( |4.47| ). We 
also make repeated use of the set of identities below which hold for momenta which satisfy 
Pi + P2 = ^3 + Pa. 

nin2{l + ^3)(1 + ^4) = (1 + ^i)(l + ^2)^3^4 

n^n2 + + ni{n's, - = (1 + ni)(l + n2)n3/(l + ^4) (4.76) 
ni(l + rii)ri3(l + n3)(iV2 - N^) = (1 + ni)(l + n2)nsni{Ns - N^) . 

Finally, since K is the momentum for an external leg we take the on shell piece: 

rk - (4.77) 

We obtain: 

N,n,{l+nk)Miil'^{K) 

= -nfc(l + ^k)NkIij^^ h — j dPdR{l + np){l+ np>)nrnkPpPppr 

NpMt-{P) , ir^i,(P)r^™(p) ^.ITUP)T%{KI 



ImSp 2 ImSp ImS^ 2 ImSp ImSfc 



(4.78) 



We introduce the symmetric notation: P := Pi ; P' := P2 ; R .= P3 and rewrite the equation 
above after symmetrizing on the integration variables. Rearranging we obtain. 
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T (K) 



^ul I dPidP2dPs{27T)^6\P, + P2-P3-K){l+n,){l + n2)nsnkPiP2P3 
,N,,Mil^{Ps) ^ N,M'i"'iK) iV,,M^^(Pi) N,,M%['^{P,) 



ImSp3 ImSjfc ImSpj ImSpj 

+ -{(iVi + iV2)S2gi) £Mgi) _ (AT, + iV3)^|4^^^^ 

+ (iVs - N^-f^^-f^ + [N, - N,)'f^'f^}] (4.79) 

Note that once again we have obtained an integral equation that is decoupled: it only involves 
Mri and rR2- With ri^2 determined by the integral equation ( [4.56| ), Equation ( [4 .791 ) can 



be solved to obtain M/j2- From the diagrams in Figs. [4,5] we see that the solutions to the 
integral equation will contain contributions of the form shown in Fig. [6]. Finally, by using 
(lOSi) , (^3TD and (PI) and comparing ( ^ISOD and (|479| ) with {^J^ and (|3|) we 



see that calculating the quadratic shear viscous response using transport theory describing 
effective thermal excitations and keeping terms that are quadratic in the gradient of the 
four-velocity field in the expansion of the Boltzmann equation, is equivalent to calculating 
the same response coefficient from quantum field theory at finite temperature using the 
next-to-linear response Kubo formula with a vertex given by a specific integral equation. 
This integral equation shows that the complete set of diagrams that need to be resummed 
includes the standard ladder graphs, and an additional set of extended ladder graphs. Some 
of the diagrams that contribute to the viscosity are shown in Fig. [7]. 




Fig. [6] : Some of the ladder and extended-ladder diagrams that contribute to the 

four-point vertex. 
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Fig. [7]: Some of the ladder and extended ladder diagrams that contribute to 
quadratic shear viscous response. 



V. CONCLUSIONS 

We have studied nonlinear response using two different methods. The first method uses 
standard transport theory. We start from a local equilibrium form for the distribution func- 
tion and perform a gradient expansion. We calculate the quadratic shear viscous response 
coefficient by expanding the Boltzmann equation and obtaining a hierarchy of equations 
that can be solved consistently. The second technique uses response theory. We work with 
a perturbative Hamiltonian that is linear in the gradient of the four-velocity field and study 
quadratic response. We generalize the Kubo formula for linear response and obtain an 
expression that allows us to calculate quadratic shear viscous response from the retarded 
three-point green function of the viscous shear stress tensor. The transport theory calcu- 
lation involves the use of the Boltzmann equation which is itself obtained from some more 
fundamental theory. The response theory calculation uses the well known methods of equi- 
librium finite temperature quantum field theory and is, in this sense, more fundamental. 
However, the response theory calculation is complicated by the need to resum infinite sets 
of diagrams at finite temperature. 

At leading order, it is well known that a correct calculation of the linear response coef- 
ficient involves the resummation of ladder graphs. Beyond leading order in response theory 
it is difficult even to identify which diagrams need to be resummed. We have identified 
precisely which diagrams need to be resummed by studying the connection between the 
transport theory calculation and the response theory calculation. We have shown that cal- 
culating the quadratic shear viscous response coefficient using transport theory by keeping 
terms that are quadratic in the gradient of the four-velocity field in the expansion of the 
Boltzmann equation, is equivalent to calculating the same response coefficient from quan- 
tum field theory at finite temperature using the next-to-linear response Kubo formula with 
a vertex given by a specific integral equation. This integral equation shows that the com- 
plete set of diagrams that need to be resummed includes the standard ladder graphs, and 
an additional set of extended ladder graphs. 

There are several directions for future work. It has been shown that the Boltzmann equa- 
tion can be derived from the Kadanoff-Baym equations by using a gradient expansion and 



keeping only linear terms [^]. The connection between this result and the work discussed in 
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this paper can probably be understood by studying the dual roles of the gradient expansion 
and the quasiparticle approximation. In addition, it would be interesting to generalize this 
work to gauge field theories. 
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